Phase Diagrams of the Harper Map and the Golden Staircase 
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We present phase diagrams of the Harper map, which is equivalent to the problem of Bloch 
electrons in a uniform magnetic field (Azbel-Hofstadter model). We consider the cases where the 
magnetic flux u> assumes either the continued fraction approximations towards the golden mean or 
the golden mean itself. The phase diagrams for rational values of u show a finite number of Arnol'd 
tongues of localized electronic states with rational winding numbers and regions of extended phases 
in between them. For the particular case oi LO = we find an infinite number of Arnol'd 



tongues of localized phases with extended phases in between. In this case, the study of the winding 
number gives rise to a Golden Staircase, where the plateaux represent localized phases with winding 
numbers equal to sums of powers of the golden mean. We also present evidence of the existence of 
an infinite number of strange nonchaotic attractors for e = 1 in points analogous to critical points 
in the pressure-temperature phase diagram of the water. [OUTP-00-07S] 

Pacs: 03.65.Sq,05.45.-a,71.30.-|-h Keywords: Azbcl-Hofstadter, Harper equation, Harper map, phase transi- 
tions, strange nonchaotic attractors, fractal, golden staircase 



The Schrodinger equation for one electron in a two- 
dimension periodic lattice on a uniform magnetic field, 
in the tight binding approximation and Landau gauge, is 
known as the Azbel-Hofstadter problem The dis- 

cretized version of this model is the Harper equation H 
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where ipn is the discretized wave function, e is the lattice 
potential strength and uj is the magnetic flux per unit 
cell in units of the quantum flux (for reviews see 
Under the transformation of variables x„ = V'n-i/V'n we 
obtain the Harper map Q 
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where (pn — {(f>o + ojn) mod 1. We take cj) £ [0, 1[ due to 
the cosine periodicity. 

The Lyapunov exponent related to the (j) dynamics is 
and the one concerning the variable x is given by 
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where 
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Prasad et al. have studied in detail the Localized- 
Extended transition for uj — u;* — (-v/S — l)/2 (golden 
mean), e = 1 and E = 0,±2.597... (center and band 
edges of the spectrum) . They found interesting connec- 
tions between strange nonchaotic attractors (SNA) ||] 
and localized-extended transition in the electronic wave 
function. We extend their analysis for all values of E and 
other typical rational values of uj. 



Our study is based on the analysis of the Lyapunov 
exponent A as a function of the parameters e, E and lu. 
Aubry and Andre ^ proved that the Lyapunov exponent 
A is proportional to the localization length 7: 
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For the Harper map A < 0. If A = 0, the phase is Ex- 
tended (E), but for A < the phase is Localized (L). 

All the characteristics discussed here are numerical re- 
sults. We start our analysis of the Harper map by dis- 
cussing the phase diagrams e x E for successive approx- 
imations ojm = Fm-i/Fm to the golden mean, where 
Fo = i^i = 1 and Fm — Fm-i + -Fm-2 (Fibonacci se- 
quence). In this paper we show the results for = 5/8 
and 8/13, which display the essential features of the 
phase diagrams for other orders of approximation for the 
golden mean. Then, we exhibit the phase diagram for 

LU = OJ . 

The e X E phase diagrams are symmetric by reflection 
in relation to the axes E — and e = 0. Then, we will 
refer only to their first quadrant. For a given a;,„, there 
are Fm-i{Fm-i — 1) connected regions of L(E) phases. 
The particular stability limits of each phase depend only 
on the value of the initial condition for the variable 4>, 
namely 0o. This behavior gives rise to co-stability re- 
gions, where more than one configuration is stable. The 
stability limits for different (po are invariant under the 
transformation 0o — > 0o + l/Fm. 

Fig. 1^ shows the phase diagram for w = 5/8. We ob- 
serve a finite number of Arnol'd tongues, a new feature in 
phase diagrams, as far as we know. The white regions (L 
phases) are labelled by winding numbers jl^ fl = p/Fm, 
where p = 0,1,2, Fm-2- In the lighter shaded re- 
gions we can find either L or E configurations, depending 
on the value of (f>o- In the darkest shaded regions there 
are E phases. The dashed curves represent the stability 
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limits of the Localized phases for (j)Q = 0, whereas the 
continuous lines represent these limits for 00 = 1/16. In 
the E = axis there are two cases. For even Fm, there 
is an E phase for < e < ei < 1; for ei < e < £2 there is 
an E /L co-stability region. For e > 62 there is a L phase. 
For odd F„i, there is an E phase if < e < 1 and a L 
phase for e > 1. In the limit m — > 00, ei and £2 — > l,and 
the behavior for odd and even Fm become the same [Q. 
3 




FIG. 1. Phase diagram for u — 5/8. The labels represent 
the winding numbers fl of Localized phases. In the white 
regions there are Localized phases, in lighter shaded regions 
there are co-stability of L and E configurations and in the 
darkest shaded regions there are Extended phases. 

The attractors of the Harper map corresponding to L 
phases, for a given Wm, are Fm — cycles, i.e., cycles with 
period Fm- In the particular case of Fig. |l|, they are 
8-cycles. The configurations corresponding to Extended 
phases are unidimensional. For uj — 5/8, typical E and 
L configurations are represented in Fig. 0. 

(a) "^(fe) 



FIG. 2. E and L configurations for lj — 5/8 and e = 0.5. 
The E-L transition takes place at _E = 0.370773... (a) Ex- 
tended phase, E = 0.34; the thick lines represent the config- 
uration and the thin lines the probability distribution of the 
iterations, (b) Localized phase, E — 0.39; this attractor has 
winding number SI = 3/8. 

In general, for a given tUm = Fm-i/Fm, the configu- 
ration of an Extended phase is made up by Fm vertical 
lines in the plane tanh(a;) x cj). The E-L transition oc- 
curs when the most visited points in the unidimensional 
E phase turn to be a Fm — cycle of the Localized phase, 
as the energy and/or e varies. 

The Localized phases within the white regions in Fig. |l| 
are actually co-stability regions. There are many stable 
Fm — cycles within each Arnol'd tongue, all of them with 



the same winding number. Fig. |^ represents the Lya- 
punov exponent A as a function of the energy, for three 
different values of 0o- We see in Fig. ^ the co-stability 
of an Extended configuration with two Localized attrac- 
tors. In Fig. ||b we see three stable attractors, all with 
winding number il — 3/8. 
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FIG. 3. The Lyapunov exponent A vs. E for lo = 5/8, 
e — 0.95 and for (j}o — (continuous line), (f)o = 1/32 (dashed 
line), 00 = 1/16 (long dashed line), (a) shows co-stability 
between Extended and Localized configurations; (b) shows 
co-stability between different Localized attractors. 

The winding number changes with the energy. Fig. ^ 
shows the variation of f2 as a function of E for uj — 5/8. 
We see plateaux with winding numbers 3/8, 2/8, 1/8 and 
0/8, which represent Localized phases. 




FIG. 4. Typical staircase for rational values of uj, represent- 
ing the winding number as a function of the energy. Here, 
u) — 5/8, e = 0.8 and (j>o = 5/8. There are four plateaux 
corresponding to the winding numbers 3/8, 2/8, 1/8 and 0/8. 
The winding number for _E = is also 3/8. These plateaux 
are connected through continuous lines, representing the Ex- 
tended configurations. 

The next approximant of lu towards uj* is 8/13. The 
number of Arnol'd tongues, as well as connected Ex- 
tended regions have increased. The stability limits for 
00 = are plotted in Fig. ||, but the full phase diagram 
must take into account all initial conditions. Neverthe- 
less, the structure of tongues with pure Localized phases 
surrounded by L-E co-stability regions is the same as in 
Fig. 1^. The winding numbers are labelling the Localized 
regions, made up by 13-cycles. 

The number of phase transition lines increases as we 
approach uj* , becoming infinite in the limit uj ^ uj* . The 
phase diagram of the Harper map for oj = uj* is in Fig. 
The labels within the Arnol'd tongues are the winding 
numbers of the Localized phases. 
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FIG. 5. Stability limits of the Localized and Extended 
phases for lj = 8/13, (^o = 0. 

The L-E transitions do not depend on the initial con- 
dition, then there are no co-stability regions. There is 
only one attractor for a given point inside a L phase, as 
can be checked by calculating the Lyapunov exponent. 
This calculation also shows the fractal character of this 
phase diagram pT| ]. 
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FIG. 6. Phase diagram of the Harper map for uj = lo* . 
The labels of the L phases are their winding numbers. The 
L-E transitions, represented as the boundaries of the L phases, 
end in SNAs for e = 1. Unlike the phase diagrams for rational 
uj, there are no co-stability regions and within each Arnol'd 
tongue there is only one attractor. The tongues to the left 
of the largest has variable winding number. Here we define 

From now on a winding number given by a sum of pow- 
ers of UJ* , say LJ*^'^ + ijj*^'^ + uj*^^ , wiU be represented by 
(fci, A:2, fca). 

In the case of Localized phases, the winding num- 
bers for Lj — oj* are limit values of those for = 
Fm-i/ Frn- In general, the phase with winding number 
Q — w* — {k) {k integer) is the limit of the phases with 
^ =^ Frn-k/ Fm, whcu m — > oo. Morcovcr, the L-E phase 
transitions end at the line e = 1. Their endpoints have 
A = and are strange nonchaotic attractors like the one 
in {e, E) = (1,0), thoroughly analyzed in reference 0. 



Therefore, there is evidence of the occurrence of an infi- 
nite number of SNAs in this phase diagram. One of them 
is shown in Fig. |^. 




FIG. 7. Phase portrait of the SNA for u = lj* and e — 1 
a,t E = 2.3447775 . . ., the endpoint of the left boundary of the 
Localized phase with winding number uj*'^ in Fig. ^ 

For e > 1 all Lyapunov exponents become negative al- 
though they maintain their bumpy structure (see Fig. |^). 
This means that the transition between different L re- 
gions is smooth without crossing an E region. This re- 
sembles the liquid- vapor phase transition of the water. 
Above the critical point the change between liquid and 
vapor is smooth. 
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FIG. 8. Lyapunov exponent as a function of the energy, for 

e — 1.1 and to = to*. The L phases occupy all the line, since 

A < 0. 

The typical E and L attractors are also different for 
UJ = Lo* . In this case, the L-E transition is between a 
unidimensional and a two-dimensional attractor, as dis- 
played in Fig. 1^. The change in the energy and/or e can 
cause a spread in the unidimensional L attractor towards 
a two-dimensional E phase, in the plane tanh(a;) x 0. 
(a) (&) 




FIG. 9. Typical Localized and Extended configurations for 
= oj* and e = 0.2. (a) E = 0.938. (b) E = 0.94. The locus 
of the L-E transition is E ^ 0.9391 .... 

One of the most original features of the phase diagram 
for UJ = UJ* is the irrationality of the winding numbers of 
the Localized phases within the Arnol'd tongues. They are 
sums of powers of the golden mean. This is the first phase 



3 



diagram to show it, as far as we know. Starting from the 
left to the right, the largest phases have winding numbers 
such that n = oj*-' , where j = 2, 3, 4, . . . , oo. The plot of 
n vs. energy for e — 0.8 shows a number of plateaux, cor- 
responding to L phases. We name it Golden Staircase 
(see Fig. ^0|) and its structure is completely different from 
the Farey tree observed in many modulated and dynam- 
ical systems 
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FIG. 11. Partial hierarchy of the L phases. 

We conjecture that the lock-in of phases with irrational 
winding numbers should be a common feature for other 
irrational values of uj for the Harper map Besides, 
the structure of the diagram f2 vs. energy should resem- 
ble the one of the Golden Staircase. 

In this letter we have presented phase diagrams for 
the Harper map forw = 5/8, 8/13 and lj* . We described 
in detail the occurrence of finite or infinite number of 
Arnol'd tongues, strange nonchaotic attractors and the 
L-E transition in each case. We also reported the occur- 
rence of co-stability regions and the peculiar structure of 
the staircases representing the winding number. 

The study of the density of states with irrational wind- 
ing numbers, following Thoulcss will be subject of 
forthcoming publication. 
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FIG. 10. (a) Winding number as a function of energy, for 
o! = CJ* and e — 0.8. The plateaux correspond to L phases and 
their winding numbers are displayed. They are sums of pow- 
ers of the golden mean. Between two plateaux there are other 
sets of plateaux, displaying a fractal structure, (b) Detail of 
the golden staircase for e = 1. 

There is an order beneath the irrational labels of the 
L phases. Firstly, the main sequence of integer powers 
of the golden mean (k) , for k — 2,3,4, .. . occurs in a 
similar way between two elements of the same sequence. 
For example, take fc = 3 and fc = 4. These labels can be 
rewritten as (3) = (4, 5) and (4) = (4, cxd). Between them 
there are the winding numbers (4, j), for j = 6, 7, 8, . . .. 
Secondly, given two neighboring tongues labelled by, say, 
(fc) and (fc -I- 1) there is a rule to find out which are the 
two largest tongues between them. For example take 
the two largest tongues between (3) and (4): (4,6) and 
(4, 7). They correspond to the lowest powers of cu* to be 
summed up to (4). It is now easy to go one step down 
in the hierarchy; take the powers (4, 6) and (4, 7), which 
can be rewritten respectively, as (4,7,8) and (4, 7, oo). 
Between them we find the winding numbers (4, 7, j), for 
j = 9,10,11,..., where (4,7,9) and (4,7,10) are con- 
nected to the largest tongues. The other steps can be 
deducted following the same reasoning. 
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